Analytical treatment of the dHvA frequency combinations due to chemical potential 
oscillations in an idealized two-band Fermi liquid 
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de Haas-van Alphen oscillation spectrum is studied for an idealized two-dimensional Fermi liq- 
uid with two parabolic bands in the case of canonical (fixed number of quasiparticles) and grand 
canonical (fixed chemical potential) ensembles. As already reported in the literature, oscillations 
of the chemical potential in magnetic field yield frequency combinations that are forbidden in the 
framework of the semiclassical theory. Exact analytical calculation of the Fourier components is 
derived at zero temperature and an asymptotic expansion is given for the high temperature and 
low magnetic field range. A good agreement is obtained between analytical formulae and numerical 
computations. 
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I. INTRODUCTION 

An ongoing question regarding magnetic oscillations in 
two-dimensional (2D) multiband metals deals with fre- 
quency combinations that are observed in de Haas-van 
Alphen (dHvA) and Shubnikov-de Haas (SdH) spectra 
even though they are forbidden within the semiclassi- 
cal theory of magnetic breakthrough (MB) by Falicov 
and Stachowiak [l|. In that respect, the organic charge 
transfer salts of the n- phase constitute experimental re- 
alizations of the linear chain of coupled orbits introduced 
by Pippard 2] in the early sixties. The Fermi surface 
(FS) of these compounds is composed of one closed tube 
(yielding the a-orbit in magnetic field) and two quasi 
one-dimensional sheets. This FS originates from the hy- 
bridization of one 2D tube, yielding at high magnetic 
field the MB-induced /3-orbit with a cross section equal to 
the first Brillouin zone area. Although these compounds, 
in particular /^-(ET)2Cu(SCN)2 (where ET stands for 
bisethylenedithia-tetrathiofulvalene), have been exten- 
sively studied, the physical origin of some of the observed 
combination frequencies % 0, p remains unclear. Actu- 
ally, besides MB-induced closed orbits which account for 
some combination frequencies such as F/3+Q, or F2/3-Q;, os- 
cillation of the chemical potential J., ^] and MB-induced 
field-dependent broadening of the Landau levels (LL) 
[3, have been invoked in order to interpret the " for- 
bidden" combination frequencies such as F/3_ck or F/3_2a 
observed in dHvA data. When they are observed in SdH 
spectra of n- phase salts 3^, ^ and of other organic met- 
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als y\ or layered semiconductor structures [T^ with sim- 
ilar FS, these latter combination frequencies can be ac- 
counted for by quantum interference (QI). Nevertheless, 
as also pointed out in the case of organic metals whose 
FS constitutes a network of coupled electron and hole or- 
bits [ll|, MB-induced LL broadening and (or) chemical 
potential oscillation as well as dimensionality certainly 
play a role in the SdH oscillatory spectra. 

The respective contribution of these two phenomena to 
the oscillatory spectra, that certainly depend on param- 
eters such as the temperature, magnetic field, effective 
masses, scattering rates, etc., needs to be quantitatively 
determined. In particular, it is of primary importance to 
quantify the temperature dependence of the amplitude of 
the various oscillations series observed in order to check if 
there is a temperature and magnetic field range in which 
the Lifshits-Kosevich (LK) formula 12] still satisfacto- 
rily applies and if addition rules such as that reported by 
Falicov and Stachowiak 1] for three-dimensional FS can 
be derived. 

In the following, a two band Fermi liquid with no 
MB and no electron reservoir (although it can have 
a significant influence on both the amplitude and the 
wave form |13|]) is considered |1J]. In such a system, 
only the chemical potential oscillation have a significant 
contribution to the dHvA spectrum. LK and Falicov- 
Stachowiak models, which assume a constant chemical 
potential [grand canonical ensemble (GCE)], predict that 
the dHvA Fourier spectrum contains only the frequencies 
linked to each orbit and their harmonics. In this case, the 
two bands are physically independent. However, for a 
constant number of electrons [canonical ensemble (CE)], 
the system lowers the total energy at zero temperature by 
sharing the electrons between the lowest LL of the two 
spectra as the magnetic field varies. In one sense, the 
filling of the lowest LL implies an effective interaction 
between the two bands, even though MB is not required. 
This leads, as already demonstrated, to a mixing of the 



fundamental frequencies |a,[i3 ^^^ ^^^ problem becomes 
less trivial. 



II. MODEL PARAMETERS 

The energy dispersion of the two bands a = 0, 1 is 
characterized by band bottoms Aq, and effective masses 
m* . In a magnetic field B^ The LL of each band are 
written as: 



ec.(n) = A, + (h'/2^ml){eB/h){n + 1/2). (1) 

In zero field, the electron densities Uoc of each band 
satisfy the conservation equation ^^ tIq, = ng, where ng 
is the total electron density. The area of the FS piece 
linked to the a-band, 5*^ = 47r^/nQ,, is given by S^ = 
27Tm^/h {Ep — Ack) or, equivalently, Sa = eFa/h where 
Fa is the Fourier frequency linked to the a orbit. 

In the following, the energies are expressed in units 
of 27Th /ttiq and we use the dimensionless variables x = 
Ue/b^ where b = eB/h and c^ = mg/m* (i. e. cq = 1 
and ci = rriQ/ml). In particular we will consider the 
case where ci is the ratio of 2 irreducible integers p/q. 
In the new units, the energies are expressed as ea{n) = 
^a + Cab{n + 1/2). It is also convenient to define the 
dimensionless fundamental frequencies fa = Sa/^Tr'^rie, 
in which case /o + /i = 1. Also defining /^ = (Aq — 
Ai)/cine, we can express the fundamental frequencies 
as a function of the effective mass ratio ci and fij: 
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Assuming a perfect crystal, the amplitude of the 
Fourier components is only driven by the thermal 
fluctuations through the damping factor Ra,p{t) — 
pAa/sinh(pAQ,), where Aq, = 27r^t/cQ,6 and t is the tem- 
perature in unit of 27rh^ /kBrriQ. For real crystals with fi- 
nite scattering relaxation time, the product Ra,p{T)R^p 
should be considered instead (the Dingle damping factor 
is given by R^^ = ex.p{—27r'^ptD/cab)^ where t^ is the 
Dingle temperature). The relations between the tem- 
perature, magnetic field and the reduced parameters are 
given by: 
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where ao is the unit cell length. In the case of the n- 
(ET)2Cu(SCN)2 compound (ag = 1.11 x l{)-^^m^,ml = 
3.2, ml = 7), the ratio t/b is equal to 2.38xT(K)/B(T). 

The chemical potential is solution of an implicit equa- 
tion which cannot be solved analytically, except at zero 



temperature (see Section IV). Nevertheless, assuming t/b 
is high enough for the chemical potential oscillations are 
weak, an asymptotic expansion can be derived at the low- 
est order. The deduced expressions allow for the calcu- 
lation of the Fourier components of the dHvA spectrum 
as reported in Section III. In the LK and FS theories 
their amplitude, which is governed by the thermal damp- 
ing factor, exponentially decreases as the t/b ratio, the 
harmonics' order and the effective masses increase. Al- 
though such feature is also observed in most cases in the 
CE, it is demonstrated hereafter that the factor in front 
of the thermal damping factor is generally not the same 
as in the LK theory. This may lead to a non monotonous 
temperature dependence of the amplitudes. E. g., the 
Fourier amplitude of the harmonics may even vanish at 
a finite temperature depending on the effective masses 
and the harmonic order. This behavior is related to the 
complex oscillations of the chemical potential that im- 
plie fluctuations of the area, hence frequencies, of the 
electronic trajectories defining the FS. More precisely, 
these frequencies are proportional to the chemical poten- 
tial /i — Aq, and therefore oscillate. It is then necessary 
to define the problem in terms of zero field frequencies, 
the oscillations of the chemical potential modifying only 
the amplitudes. 

The analytical formulae deduced in sections III and IV 
are compared to the Fourier spectra deduced from nu- 
merical computations of the chemical potential both at 
zero and finite t/b. Discrete Fourier transforms are cal- 
culated with a Blackman window in a given field range 
from bmin to bmax • The absolute value of the Fourier am- 
plitude {A) at a given frequency for a mean field value 



2/(1/6, 



^/bmax) is determined from the 



amplitude of the discrete Fourier transform (Acaic) as 



A = 4A,alc/0.84{l/br, 



^/braax)' As reported here- 



after, an excellent agreement between analytical expres- 
sions and numerical computations is obtained at zero 
temperature while the asymptotic expansion is a good 
approximation in the large t/b range. 



III. FOURIER COMPONENTS OF THE 

MAGNETIZATION IN THE HIGH 

TEMPERATURE REGIME 

To evaluate the thermodynamical quantities, we be- 
gin by considering the Fourier transform of the grand 
potential Vt in the GCE for the two band system. It is 
expressed as: 
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The Fourier transform of the magnetization M = 
—dVt/dB at fixed chemical potential ji contains cosine 
functions of the arguments 27rpFa/B = 27rpfaX with p in- 
teger, and therefore the frequencies involved in the spec- 
trum are the individual frequencies Fa plus their har- 
monics. This is the LK result. In the GCE, there is 
therefore no mixing of the form kfo + Ifi with (/c, /) pos- 
itive or negative integers. In the CE, the magnetization 
is expressed as M = —dF/dB, where F = Q-\- riefi is the 
free energy. Contrary to the GCE case, it is known that 
the Fourier spectrum of the magnetization contains such 
kfo + Ifi combinations, where the fa = (/io — ^a)/ca^e 
are defined relatively to the zero field chemical potential 
/io = /J^{B = 0). The chemical potential is solution of the 
equation: 



Efi- Aa 
— =ne 



(3) 



-^EE 



^^R^^p{t) sin 27rp^——^ 
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At zero field, the last term vanishes, and we get the 
simple relation jj.{B = 0) = jhq = Aq, + CaUa- At finite 
magnetic field, the chemical potential oscillates, and is 
solution of the self-consistent equation Eq. ^. If we 
replace /i in Eq. ((SJ by this solution, the free energy is 
obtained using the relation F = Q -\- riefi. It can be 
noticed that the series in Eq. (|3|) is proportional to the 
small parameter b. We may therefore replace /a in the 
cosine functions by its zero magnetic field value /io- 
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This relation can be expressed as: 



11^ flQ 






(4) 



with M^ = Ep>ii-^yRpAt)/^P X M^^Pfp^)' Re- 
placing /i in Eq. Q by this expression, and then com- 
puting F = Q -\- riefi^ we obtain a finite temperature 
approximation for the free energy. The approximation 
Eq. ((5J shows that the dominant oscillatory behavior of 
the potential leads to frequency combinations. Indeed, 
in the oscillatory part of Eq. (|21), the arguments of the 
cosine functions are themselves oscillating: 



cos 27rp 
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Cab 
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(5) 



X TTexp(— 2z7rpK;Q,M(3(x)) > , 
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where Wa = c^^ / ^g Cq^ are mass weights. The func- 
tions e-x.-p{—2iT:pWaMi3{x)) are periodic functions of x 
with ff3 period, therefore they can be expanded in their 
Fourier series: 



exp(-2i7rp^c.M;3) = ^ 5c,,/3(p, g) exp(2z7r(7//3x). (6) 

q— — oo 

We obtain a product of several Fourier series, and ex- 
panding the product and combining the oscillating func- 
tions lead to the combinations between the fa frequen- 
cies. 

If we try to go beyond this approximation, the scheme 
breaks down. Indeed, 5 ja = ja — jiiQ is di small parameter, 
and we would like to expand the sine functions in Eq. (|3|) 
relatively to this parameter. The first order is the case 
we discussed before, and at the next order we obtain 



sin 27TP 
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Cab 
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COS {27rpfax) Sjl. 



(7) 



The last term in Eq. Q is proportional to p d/j^/b o:^ p 
which is not small. It is therefore difficult to expand 
Eq. Q) as a series of dfi. However, the first order approx- 
imation works for high t/b values (this is checked here- 
after by solving numerically the self-consistent relation 
Eq. 0), since the oscillations of the chemical potential 
are negligible in this range. The oscillating part of the 
magnetization is defined in units of he/rriQ by rriosc oc 
{ne/b'^)dF/dx or rriosc = T.f ^ce{F) sm{27rFx). The 
dominant value is: 
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^ BaAPa^ Qo) exp{2i7rqofox) 
^ Ba,i{Pa, qi) exp{2i7rqifix) 



(8) 



The various terms in Eq. ((S|) comes from the derivation 
of Ue/J^ + ^ using Eq. (0) and Eq. Q). In the large t/b 
range, the coefficients Ba^p can be estimated by noticing 
that Mp is dominated by the first term of its Fourier 
series M(^{x) c^ — i?/3^i/7rsin(27r//3x). In this case those 
coefficients are simply related to Bessel functions: 



BaAP^q) - Jq(2pWaR(3,l)' 



(9) 



The arguments of the Bessel functions are rather small 
since Rf^^i '^ 2A/3 exp(— A/3) <C 1. The amplitudes 



^ce{F) and the related masses can be extracted within 
this approximation at the lowest order, and for different 
frequency combinations. The m^^ harmonics of /a, result 
from two contributions: one coming from the product of 
Ma{x)M^{x) in Eq. 0, and the other from the solutions 
{Pa^Qa = ^rna,qf3y^a = 0) and {pp^qp = 0,qa = ±ma) 
in the second term of the right hand side of Eq. (|H|) . For 
rria = 1, we simply obtain: 
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47rne/a- exp(-27r^ — -) 

b Cab 



(10) 



which is the same amplitude as in the GCE. In- 
deed the main contribution comes from the terms 
B a,a{^^^)Ba, f3{^^0) c:^ 1, the ucxt Corrective terms are 
smaller and correspond to the third order in Ra^i- The 
amplitudes of the second harmonics AcEi'^fa) are differ- 
ent in the CE: 



AcEi2U) ^ Agce{Vo.) 



^GCE 



(2/.) 
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1 - 47r — -Woe 

Cab 

-47rne/a-exp(-47r^ — -). 
b Cab 



(11) 



Within this approximation, AcE^^fa) vanishes and 
changes sign at t/b = CQ,/(47r^i^Q,), due to the additional 
factor in Eq. (|TT|) . Additional corrective terms may mod- 
ify this value, but a more explicit expansion would be 
needed in order to get a better accuracy. 

The amplitudes for the sum and difference of /o and 
/i are computed in the same way and we obtain: 



^CE[ 
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The amplitudes for f^ + 2/^3, a ^ j3^ are 



AcEiU + Vp) ^ ^^^fe^lr^^-.i X 



( Rl3,2 , ,Wa , x-o2 \ 

V^ + ^T + ^'^^M 



(13) 



leading to: 



AcEJU + 2/g) ^ /^ + 2/g 1 + (1 + W0)Xp 

AcEiffi + 2U) ~ ff} + 2Ul + {l + Wo)Xa ' 



(14) 



The amplitudes given by Eqs. ^J to ^J are generally 
the largest that can be extracted from data analysis. The 
ratio 



Ace {fa - ffd) f(3 - fa 



AcE{fa^ff3) /a + //3' 



(15) 



is constant at high temperature indicating that these 
two oscillations have the same mass m^{fa^ff3) = m^ + 
rrip that can be extracted by plotting the logarithm of 
{b/t)'^AcE{fa^f(3) versus t/b^ which should be a straight 
line in this limit. More generally the addition rule 



m*{kfa + lfp) = \k\ml 



rrifi 



(16) 



can be derived. Nevertheless, it should be kept in mind 
that the t/b dependance of the prefactor is not simple. 
Therefore, each case has to be treated separately. E. g., 
in the high t/b range, Eq. E|can be rewritten as: 



AcE{fa^2fp) 



4ne(/c + 2//3)Ac,A/3 



X [1 + 2A;3(1 + wp)] exp[-(A« + 2Xp)]. (17) 

This leads to m*(/c + 2//3) = m* + 2m^, as predicted 
by Eq. ^J although this effective mass is derived in 
the high temperature range by plotting \og[AcE{fa + 
2fi3){b/t)^] vs. t/b. Oppositely, the asymptotic factor 
in front of the dominant exponential damping is always 
proportional to t/b in the GCE. 

The asymptotic expressions Eq. (|10|) to Eq. (|13|) (see 
the solid lines in Figs. Q and 12) are compared to the 
amplitudes of the relevant Fourier components deduced 
from the direct numerical resolution of the chemical po- 
tential given by Eq. ^ (symbols in Figs, ^andl^l). First, 
the non monotonous behavior of Ace {2 fa) predicted by 
Eq. (|TT|) can be observed in numerical data of Fig. QJ) 
although at a shifted t/b value. More generally, a good 
agreement between numerical and analytical data is ob- 
tained for t/b > 0.05 and t/b > 0.1 {t/b > 0.15 for /o+2/i 
and 2/o + /i) in the GCE and the CE, respectively. This 
is also evidenced in Fig. O in which the magnetization 
oscillations deduced from numerical resolution of Eq. ^ 
at t = 4 X 10~^ and t = 7 x 10~^ are compared to the 
analytical approximations given above, in the field range 
0.02 <b< 0.05. 



IV. ZERO TEMPERATURE FOURIER 
ANALYSIS 

At low temperature, the amplitudes differ strongly 
from the LK theory, due to the importance of the chem- 
ical potential fluctuations. At zero temperature, the 
chemical potential equation Eq. ^ can be solved exactly 
and expressions for the total energy and magnetization 
are obtained. The total energy £^ is a function of x and 
can be expressed as 
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FIG. 1: Reduced temperature (t) and magnetic field (b) de- 
pendence of the dHvA oscillations amplitude for the canonical 
ensemble. The band parameters are Ao — Ai = 1/10, c^^ = 
1, c-^ = 2 (leading to fo = 4/15 and fi = 11/15). Symbols 
are data deduced from numerical computations in the reduced 
temperature range from t = to t = 0.009. Solid lines are 
deduced from the high t/b approximation given by Eqs. (fTn|) 
to (CSl. Note that, Eq. ^ predicts that A{2fo) and A{2fi) 
vanish at t/b — Ca/{4:7T^Wa)- 



E{x) = Eo ■ 



^G(x), 



(18) 



where G{x) is an oscillatory function, and Eq is the 
energy in absence of field. The latter is given by 



Cl 



2(1- 

(Ao 



Cl) ^ 
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-(ciAo + Ai) 



2(1 + ci) 



(19) 



In the following, we compute G and its Fourier trans- 
form by using a combinatorial analysis. We give here the 
main analytical arguments to obtain the final result. If 
we suppose that Aq is large compared to Ai (the final 



FIG. 2: Same as FigQ for the grand canonical ensemble. 
Solid lines are deduced from Eqs. (|1()|) and (|11|) . 



result does not depend on this assumption), we first be- 
gin to fill the LL corresponding to the band a = 1 since 
it has a lower band bottom. If ni{x) is the number of 
LL we have to fill before reaching the first level of the 
band a = 0, then there are b{x — ni{x)) = bxy electrons 
left we have to share between the two bands. We obtain 
^i(^) = [(Ao — Ai)/cib — 1/2], where [. . .] means the 
integer part. We can relabel in this case the LL energies 
of the 1-band as ei(n) = Ai + cib{n + 1/2), n > 0, with 
Ai = Ai + cib[{Ao - Ai)/cib - 1/2] and Ai -^ Aq when 
b ^ 0. The LL corresponding to n < of the 1-band are 
already filled with bni{x) electrons and there are now 
bxy electrons for the remaining LL in the bands and 
1. If we assume that Ci = p/q is a rational with (p^q) 
positive integers, then we can truncate the total Landau 
spectrum, made of the superposition of the LL of the two 
bands, by noticing that this spectrum has a periodicity 
equal to p-\- q. We call periodicity the minimum number 
of LL after which the relative positions of the LL between 
the bands 1 and 2 repeat themselves. We can therefore 
decompose x as x = ni {x) -\-ko{p^q)^ mo + mi + r , with 
/co, mo and ttii integers, ko is the maximum number of 
blocks ofp-\-q LL completely filled, m^a is the number of 
LL in the a-band filled, and < r < 1 is the fraction of 
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FIG. 3: DHvA oscillations deduced from numerical compu- 
tations for the canonical (CE) and grand canonical (GCE) 
ensembles, at two different temperatures (t = 0.004 and t = 
0.007). The band parameters are the same as in Fig. Thin 
lines are deduced from Eqs. (|1()|) to (|13|l . 



the last LL occupied by the remaining hxr electrons (ei- 
ther in band or 1 in accordance with the constraint that 
the total energy should be minimum). We have the rela- 
tions < rriQ < p and < mi < ^, < mo + ttii < p-\- q. 
ko is simply defined by ko = [y/{p -\- q)]^ r = x — [x] and 
the sum mo + mi = [y — ko{p + q)] = M. We finally have 
to find either mo or mi, since the sum is determined, and 
the last LL occupied by the b x r electrons. The total 
energy is 



Al)/2ci (the energy corresponding to the LL of the 1- 
band filled between Ai and Ao). We obtain mo or mi 
by noticing that they are independent of r. Therefore, 
we compute £^ as a function of, e. g., ttio for r = 0, and 
minimize the energy with respect to m^Q. We obtain a 
quadratic function of ?no, and the minimum is given by 
the integer closest to the bottom of the parabolic curve: 
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(21) 



The energy can be put diS E = E^ -\- ipEcik^^ M, r)/2. 
The effective inverse field x = Ue/h (or y) can be varied 
by either changing the magnetic field h or changing the 
density of electrons ng. In the latter case, all quantity 
depending only on h is constant, such as £^5 for example. 
This is especially useful if we want to extract the peri- 
odical part of E on the top LL. The periodicity of the 
energy is directly related to the periodicity of the LL, 
and this happens every time that {p -\- q)b electrons are 
introduced in the system, hence at fixed magnetic field 
b. Following this argument, the ko dependence in Ec can 
be removed by iteration, if we express Ec{ko) as function 
of Ec{ko — 1). We obtain 



Ec{ko,M,r)=Ec{0,M,r) 



(22) 



'' -(^2-(M + r)^) + 2^--^^(ciAo-Ai). 
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We then separate the terms depending on M and r 
and the terms depending on y: 
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We can check that the sum of the first 3 terms on the 
right hand side of the equation Eq. (|23|) are equal to £^0 
[see Eq. (|T9J) ] with some additional terms: 
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The first term gives the bottom energy E^ = (A^ 
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We can therefore identify the G function in Eq. I|18(l 
as a periodic function of the variable y: 



G{y) = moiyf + Cl {[y] - mo{y)f 
+ ci5{[y] - 2mo(y)) + {y - [] 



(25) 



mo{y) + 2 + 



ci{[y] -mo{y)^-) 
ciS + ci{[y] - mo{y) + 2^ ~ ^^'^^^ ~ 2 



1 + ci 



1 + Ci 



1 + Ci 



where 5 = (Ai — Ao)/ci6, ?/ being constrained to be 
inside the interval < y < {p -\- q). Indeed, using the 
fact that mo{y -\- p -\- q) = mo{y) + p, we check that 
Eq. (|25|) satisfies G{y-\-p-\-q) = G(^). In the fohowing, we 
expand G{y) in Fourier series, and, using y = x — ni{x)^ 
express the result as a function of x, in order to obtain 
the relation Eq. (|T^ . The Fourier coefficients are defined 
by: 



G, 



G{y) 



/ 



(P+9)/2 / ly 

G{y)exp -2iiT — — ) dy 

(p+g)/2 V P + ^. 



1 "^ ~ / 

— ^ G^ exp f 2i7r 



p + q 



l= — oo 



ly 

p + q 



(26) 



Defining the frequencies co = 27rl/{p + (7), the integra- 
tion over y can be split into small parts (assuming for 
example that {p + q)/2 is integer): 

b+9)/2 .1 

G= ^ exp{-iul) dr G {I ^r)exp{-iur). (27) 

In Eq. ((23, r?io(^ + r) = mo{l)^ and G is a second 
order polynomial function in the variable r, hence the 
integration over r is straightforward since it uses the in- 
tegrals Jq exp(— ia;r)r^(ir, with /c = 0,1,2. In the limit 
of irrationnal ratio ci, as p + g ^ 00, the computation 
is simpler and we assume that the physical quantities 
should be smooth functions of ci. In this limit: 



G(cj) = / G{y)exp{-iu;y)dy 

/oo 
G{u;)exp{iijy)du;. (28) 

-co 

In Eq. (|27|) the sum is now from —00 to 00 and we have 
to evaluate the following distributions 

CXD CXD 00 

2_] exp(— zo;/), 2_] exp(— zo;/)/, V^ exp(— zo;/)/^, 

l= — oo l= — oo l= — oo 

00 00 

^ exp(— za;/)7Tio(05 /. exp(— za;/)7Tio(05 

Z= — CO Z 

CO 

2^ exp {—iu I) mo {1)1 



l= — oo 



in order to obtain the final amplitudes G{lj). The sim- 
plest series is 



Y^ exp {-icol) = 271 Y^ S{u; - 2tiI) = S{uj). (29) 



I— —00 



I— —00 



The other series Xl^^_^ exp{—iul)l^ with /c = 1, 2 are 
derivatives of 5(0;), and include therefore derivatives of 
S functions. It appears that these derivatives cancel each 
other in the final result. The remaining distributions can 
also be written in term of S{uj) function and its deriva- 
tives: 



ciS 



2_\ exp{—iLjl)mo{l)(x- S{ij) 



(30) 



E 



{-ly exp{2i7rlciS/{l + ci)) ^ f 2itIci 



2M 



l + cW ' 



and 



Y exp {-iujl) mo{lf ^ f Y^ + 7Y 



2r2 



1 cfS' 



(1+Cl)^ 



S{u) 



^^ {-ly exp{2i7rlciS/{l + ci)) 



1 



2i7rl (1 + Ci) 






1+Ci 



Finally, the series Xlz^-oo e^P(~^^O^o(0^ depends 
only on the derivatives of S. The details for obtaining 
Eq. (pnj) and Eq. (pTTj) are given in appendix ^ G(cc;) 
can be expressed as a sum of delta functions at frequen- 
cies corresponding to the arguments of the function S. 
After some tedious computations, we arrive at a Fourier 
spectrum containing frequency combinations /o/o + ^i/i- 
In order to obtain G{x) as a Fourier series, we have to 
replace yhyx — no (x) in the cosine functions and express 
the result as a function of x only. Since no{x) is integer, 
we always have cos{27rly) = cos(27r/x). It is also useful 
to notice that y -\- S = x{l — fb). Finally we obtain 



G{x) 



1 1+4:01+4 



12 1 + ci 



(32) 



1 = 1 



(-1)^(1 + ci)^ / 27rki 

' 1 — cos z I COS [27tIjox) 



l= — oo 



, oW (-1)''-' {\ + Cif 

X ( 1 - COS ^^(^'-Oci \ ^^g (27r(r/o + l!x)x) . 

1 + Ci / 



Eq. ((5^ is the main result for the two-band prob- 
lem at zero temperature. It contains the exact ampli- 
tudes of the oscillating part of the energy. The oscil- 
lating part of the magnetization is simply defined by 
ruosc = neG'{x)/2 = Y.^, Ace{F) sm{27rFx). We no- 
tice that G{x) depends only on ci and the fundamen- 
tal frequencies fa- In the double sum [Eq. ((S^ ]. it may 
happen that some terms diverge because / + Vci van- 
ishes in the denominator. But in this case, the quantity 
(/ — V)ci/{1 + ci) = / is integer and the interference term 
1 — cos(27r/) vanishes at the same time. It is easy to check 
that the quantity (1 — cos(e))/e^ is finite when e tends 
to 0, and the divergences do not occur. The values of 
these singular terms are given by this regularization. In 
the general case, we can give the expression of the ampli- 
tudes AcE{kF) for the k^^ harmonics of F = /q/o + hfi^ 
where /i is a positive integer and Iq can be negative. If 
/o//i is equal to an irreducible ratio of 2 integers pi/po^ 
then we can find many combinations I' fo-\-lfi in the dou- 
ble sum of Eq. (|S^ that are equal to kF. This occurs for 
I = kli -\- mpi and I' = klo — mpo. We obtain 



AcE{k{fo^fi)) 



Cl 



Ci 7rk 



(36) 



Other interesting cases are harmonics of /i and /q. The 
amplitudes for kfi are obtained by considering Eq. (|32j) 
with I = k -\- mpi and V = —mpo: 



AcEikfi) = Y. 



m>0 



{k + m{po + pi)Y{k + m{pi - cipo)y 



1 — cos 



27r(/c + m(po+Pi))ci\ kf_ 



1+Ci 



27r3 



h (37) 



In the irrational case, only the first term m = con- 
tributes: 



Acsikfi) = 



(_l)fc+i(l + ci)Vi 
27r3fc3 



1 — cos 



27rfcci 

1+Ci 



.(38) 



AcEikF) = 



kF{l + ci] 



27r3 



E 

m>0 



f_l\k{lo-li)-m{po-\-pi)-\-l 

{k{li -h Clio) + m{pi - cipo))^ 



1 — cos 



27r{k{lo-li)-m{po-\-pi))ci 
1+ci 

{k{lo -h) -m{po ^Pi)y 



(33) 



In the case where /0//1 is irrational, only the term 
m = has to be taken into account: 



AcEikF) 






(34) 



[1 



COS 



27rk{lo — /i)ci- 
1 + ci 



The amplitudes Eq. (|33j) holds for any frequency ex- 
cept for the special case k{fo -\- fi). Indeed AcE{k{fo + 
/i)) is obtained with the condition po ^ p\. 



For the harmonics /c/o, when /0//1 is rational, we take 
/ = mp\ and V = k — mpo^ with the condition k 7^ m{po-\- 
pi), and we obtain 



AcEikfo) = Yl 



m>0 



{k - m{po +pi))^(A:ci + m(pi - cipo))^ 



,_,^,Mk-m{po^p,))cAkfo^ (39) 

1 + Cl J 27T^ ^ ^ 



If /0//1 is not a rational, only the first term m = 
contributes as before: 



AcE{kfo) = o 0.0 2 1 - cos — — - .(40) 



27r^k^cl 



1 + ci 



Usually the terms m > 1 are only small corrections due 
to the I/tti^ behavior. Finally, the "forbidden" frequency 
amplitude corresponding to /i — /o is given by taking 
/ = 1 + mpi V = -1- mpo in Eq. ((23) 



AcE{k{fo^fi)) 
k 



Cl 



E 



1 + Cl tt/c 

(l+Ci)^(-l)^(^o+^i) + ^ 



27r3(po +Pi)^ ;^^ m2(A:(l + ci) + m(pi - cipo))^ 



1 



27rm{po +pi)ci 
1 + ci 



(35) 



In the case where /0//1 is irrational, only the first term 
on the right hand side of Eq. (p^ has to be taken into 
account, and we obtain a simple result: the k dependence 
is the same as it would be in GCE case for a single band 
with fundamental frequency /o + /i (with the exception 
of a (—1)^ term that appear in the GCE case): 



AcE{fi-fo)=Yl 



(1+cl)^ 



(2 + m(po+Pi))^ 



X 1 — cos 



m>0 

f_l\'m{po-\-pi)-\-l 

y^ ^ ^ 

(1 - Cl +m(pi - cipo))^ 
27r(2 + m{po + Pi))ci \ fi - fo 



1 + Ci y 27r3 * 

For the irrational case, we just have 



AceUi - fo) 



{l+c,fh-h 
(l-ci)2 (27r)3 



1 — COS ■ 



47rci 

1 + ci 



(41) 



(42) 



These results imply that the jumps of the magnetiza- 
tion at zero temperature only come from the contribution 
of the /o + /i harmonics in the case when /0//1 is irra- 
tional (or P07 Pi infinite). Indeed, a jump in magnetiza- 
tion means that the absolute value of the k^^ harmonics 
in the free energy or grand potential Fourier spectrum 
decreases for large k like 1/k'^ (respectively 1/k for mag- 
netization). Only in Eq. (p^ and Eq. (|36|) this condition 
is satisfied. In other cases [Eq. (|33j) ]. the decreasing is 
only like l/k"^ (1/A:^), as for the frequencies correspond- 
ing to kfo [Eq. (gOl] or kfi [Eq. (|38|l]. 

Fig. El compare the Fourier analysis of the dHvA spec- 
trum deduced from numerical computation of the chem- 
ical potential [see Eq. ^] to the predictions of Eqs. 
([HIj) . ((36j) and ^^ to (gOJ (see horizontal dashes in Fig. 
0J. An excellent agreement between numerical compu- 
tations and analytical formulae is obtained even in the 
case where a given frequency results from several combi- 
nations. This is e. g. the case of the Fourier peak at F 
= 8/15 in Fig. ^ which results not only from 2/o but 
also from e. g. 13/o — 4/i and — 9/o + 4/i. Nevertheless, 
these two other Fourier components have a small ampli- 
tude when compared to Ac£;(2/o) so that formulae such 
as Eq. (|39|) which is relevant for this example, remain a 
good approximation in that case. 
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CONCLUSION 



The Fourier components of the de Haas- van Alphen os- 
cillation spectrum have been calculated for an idealized 
two-dimensional Fermi liquid with two parabolic bands in 
the case of canonical and grand canonical ensembles. Ex- 
act analytical calculations of the Fourier components at 
zero temperature for the CE [see Eqs. (IHH), ^^ and (p^ 
to (001)] are in excellent agreement with numerical com- 
putations of the chemical potential [Eq. ((HI)], as demon- 
strated in Fig. ^ Regarding the high t/b range, asymp- 
totic formulae have been derived for the CE and GCE [see 
Eqs. (|10|) to ([T^ ]. Addition rule of the effective masses 
associated to the combination frequencies observed in the 
CE has been derived (see Eq. (|TB|) ). Nevertheless, and 
contrary to the GCE case, the t/b dependence of the 
Fourier components' pre-exponential term is non-trivial. 
In particular, the 2^^ harmonics amplitude can vanishes 
and changes sign as t/b varies. A good agreement is ob- 
tained between analytical formulae and numerical com- 
putations for t/b higher than 0.15 in the less favorable 
cases (see Figs. n]to|3|. Such a value which corresponds 
to e. g. T/B > 0.06 K/T for /^-(ET)2Cu(SCN)2, is suit- 
able for the determination of the effective masses from 
experimental data with a good accuracy in a large tem- 
perature and magnetic field range. 



FIG. 4: Fourier analysis of dHvA spectra at zero temperature 
computed in the canonical ensemble case. Horizontal dashes 
are results deduced from Eqs. ElESlEHlEniand|ini The band 
parameters are Ao — Ai = 1/10, c^^ = 1, cj"^ = 2 (leading 
to fo = 4/15 and fi = 11/15) and Ao - Ai = 2/47, c^^ = 1, 
c-i = 2670/2867 (leading to fo = 52/113 and fi = 61/113) in 
Figs. 2Jl and^J), respectively. 



APPENDIX A 

The series Xl^^-oo exp(— zcc;/)mo(0 can be evaluated by 
considering the following periodic function: 



(x) = mo{x) 



ci 



ci 



-X, 



with 



mo(x) 



1 + ci 2 



It is indeed simpler for calculations to consider mo in- 
stead of mo, both functions having the same values for 
integers mo{l) = mo{l). We have (j){x -\- p -\- q) = (j){x)^ 
and therefore the Fourier transformation is written as 



{x) = 



ci5 

1 + ci 



(Al) 



E 






From that result, we have 



10 



we obtain after further computation Eq. (|3T|) . The abso- 
lute value in Eq. ((23 is periodic with period p -\- q and 
therefore the function we can consider for example is 



E 

l= — oo 



ex.p(—iujl)mo(l) = \^ ex.p{—iujl) I 0(/) H / 1 , 

, V 1 + ci y 



the second sum on the right hand side of the previ- 
ous equation only contributes to derivatives of S. The 
sum over (j){l) can be performed if we use Eq. (|A1|) and 
commute the double sum. We obtain Eq. (|30j) : 
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(l+ci)mo(a;) -ci{x+ 2) + o ~ ^^^ 



■ (A4) 



As previously it can be Fourier transformed and the se- 
ries X]^^-cx) exp(— za;/)^(/) is expressed as a sum of func- 
tions S{(jo). We obtain 
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We use the same analysis to compute the series 
X]^_oo ^-^P(~^^0^^o(0^- Here we have to consider in- 
stead the function 



(j){x) = I rho{x) 



ci 



■ci 



(A3) 



This is also a periodic function of x with period p + <7, 
and it can be Fourier transformed. As (j){l) oc mo{l)'^, 
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